Robust charge and magnetic order under electric field/current in the multiferroic LuFe204 
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We performed elastic neutron scattering measurements on the charge- and magnetically-ordered multiferroic 
material LuFe204. An external electric field along the [001] direction with strength up to 20 kV/cm applied 
at low temperature (~ 100 K) does not affect either the charge or magnetic structure. At higher temperatures 
(~ 360 K), before the transition to three-dimensional charge-ordered state, the resistivity of the sample is low, 
and an electric current was applied instead. A reduction of the charge and magnetic peak intensities occurs when 
the sample is cooled under a constant electric current. However, after calibrating the real sample temperature 
using its own resistance-temperature curve, we show that the actual sample temperature is higher than the 
thermometer readings, and the "intensity reduction" is entirely due to internal sample heating by the applied 
current. Our results suggest that the charge and magnetic orders in LuFe204 are unaffected by the application 
of external electric field/current, and previously observed electric field/current effects can be naturally explained 
by internal sample heating. 
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LuFe204 is a new multiferroic material. Here the bulk fer- 
roelectric polarization is not due to cation displacements as in 
conventional ferroelectrics, but instead arises from the three- 
dimensional charge valence order of Fe^+ and Fe''+ ions oc- 
curring at ^ 340 K. ' The magnetic order starts at lower tem- 
perature 240 K) in the charge-ordered ferroelectric phase. ^ 
Recent reports of strong couplings between the two orders''^, 
as well as large room temperature magneto-electric response 
in this material^'"*, make LuFe204 a promising candidate for 
practical applications. 

In addition to the magneto-electric response, tremendous 
interest has been focused on studying the electric-field re- 
sponse of the magnetic structures in multiferroic systems. 
Nevertheless, there have been only a few observations of 
such effects^""' in known multiferroic systems, and all of 
them can be attributed to electric field realigning ferroelec- 
tric domains and therefore causing a macroscopic magnetic 
response. In LuFe204, the ferroelectric polarization is charge- 
valence driven, and the charge-valence order also couples 
strongly to the magnetic order. If the charge order can be 
affected or broken, it is then possible to affect the micro- 
scopic magnetic structures by an external electric field or 
current. Indeed, there have been previous reports on non- 
linear current-voltage behaviors, and eventually an electric- 
field "breakdown" of the charge order in LuFe204 There 
have also been claims of electrical control of the magnetic 
response in the same material»i^ These observations, if con- 
firmed and fully understood, would be extremely interesting 
and important for achieving mutual control of electric and 
magnetic degrees of freedoms in multiferroic systems. 

We thus performed elastic neutron scattering measurements 
on single crystals of LuFe204, studying the response of the 
charge- and magnetic-order Bragg peaks under external elec- 
tric field/current. No electric-field effect has been observed at 
low temperature ^ 100 K for a field strength up to 20 kV/cm. 
Near room temperature an electric current effect on the order- 
ing is observed, and we show that it is due to internal heating 



of the sample by the current flowing through the sample. We 
conclude that the charge and magnetic order are robust and 
not affected by the electric field/current. 
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FIG. 1: Schematic of the magnetic- (open circles) and charge-order 
(filled circles) peaks in reciprocal space. The arrows indicate the scan 
directions along [110] and [001]. 

Single crystals LuFe204 were grown using floating zone 
technique.'"* Typical crystal sizes are ~ 10 x 5 x 2 mm"^. 
Our neutron scattering measurements were performed on BT9 
triple -axis -spectrometer at the NIST Center for Neutron Re- 
search. An incident neutron energy of 14.7 meV was selected 
by a pyrographic (PG002) monochromator, with beam col- 
Hmations of 40'-40'-40'-80', and another PG002 crystal was 
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used as the analyzer. PG filters were used before the sample 
to reduce background from higher order neutrons. A sam- 
ple of 0.9 g was loaded in a closed-cycle refrigerator, where 
the thermometer is attached to the mounting base, about 3 cm 
away from the sample. LuFe204 has a hexagonal structure 
with three iron double layers in each unit cell. The ferro- 
electric polarization is directly due to the imbalance of iron 
valences in each double layer, and the net (induced) polariza- 
tion appears along the [001] direction, perpendicular to the 
hexagonal plane and the double layers. The two 10 x 5 mm^ 
(001) surfaces were painted with silver paint so that electric 
field/current can be applied along the [001] direction (2 mm 
thick) for our measurements. The single crystal sample was 
oriented so that the horizontal diffraction plane is the [HHL) 
plane, defined by the vectors [110] and [001]. (See Fig.[T]). 
The resistance-temperature (R-T) curve has been measured 
using Keitheley 2000 multimeter, and the in-situ resistance 
during the neutron scattering measurements was obtained by 
reading the voltage across the sample while keeping the cur- 
rent constant. 

The magnetic Bragg peaks in this compound occur at recip- 
rocal space positions (1/3, 1/3, L) and (2/3, 2/3, L) for both 
half-integer and integer L values '^ '^ below the magnetic or- 
dering temperature Tjv ^ 240 K; while the charge peaks only 
appear at half-integer L value a^-^''^" — . (See Fig.lT]!. 
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FIG. 2: (Color Online) Scans through (a) charge-order peak (2/3, 
2/3, 3.5) along [001] direction at different temperatures; (b) magnetic 
peak (1/3, 1/3, 0) along [001] direction; (c) satellite peak (0.306, 
0.306, 0) along [110] direction. Error bars represent square root of 
the total counts, and those in (a) and (b) are smaller than the symbols. 
Lines through data are guides to the eyes. 

We choose to monitor (2/3, 2/3, 3.5) for the charge order, 
and (1/3, 1/3, 0) for the magnetic order. Representative scans 
through these two peaks are plotted in Fig.|2|a) and (b) respec- 
tively. These peak intensities indicate the charge/magnetic or- 
der and are plotted vs. temperature in Fig. [3 a) and (b). For 
zero field, the charge order starts as two-dimensional (2D) at 
around 500 K'^, and becomes three-dimensional (3D) around 
Tco ^ 340 K^iifi. At Tn 240 K, the magnetic order oc- 
curs, as indicated by the rise of (1/3,1/3,0) intensity shown 
in Fig.[3lb). Also, a boost to the intensity at (2/3,2/3, 3.5) is 
observed [see Fig. Oa)], which is due to the additional scat- 
tering intensity at this wave-vector coming from magnetic or- 
dering. The intensity at (2/3,2/3,3.5) now (for T < Tm) has 
contributions from both the charge and magnetic orders. With 
further cooling, the intensity increases until reaching another 



temperature Tl ^ 180 K. Here another phase transition oc- 
curs, similar to that observed in Ref. Note that this second 
phase transition is strongly sample dependent, and has been 
shown to be missing for some s ample Si2iiii2 As demonstrated 
by previous studies, the magnetic properties of LuFe04+5 are 
sensitive to oxygen stoichiometryrS and different temperature 
dependence of intensity below Tl ~ 1 80 K is likely due to dif- 
ferent oxygen content in different samples. The intensities for 
both peaks drop at Tl but become almost flat below ~ 150 K 
[Fig-EIa) and (b)]. This is also demonstrated in Fig.|2] where 
at 100 K both the charge and magnetic peak intensities are 
lower than those at 180 K. Additionally, at 100 K there is a 
strong 2D diffuse type magnetic scattering'^ which shows up 
as higher "background" in the L-scans [Fig.|2jb)]. In addition, 
intensity starts to appear at satellite positions below Tl, and in 
Fig.|2c) we plot scans through a satellite peak around (0.306, 
0.306, 0), whose intensity dependence on the temperature is 
plotted in the inset of Fig.|3jb). 

These results suggest that the low-temperature (mag- 
netic/charge) structures of LuFe204 are quite complicated and 
sometimes sample dependent. Nevertheless, our goal is to 
search for possible electric-field effect on the magnetic and 
charger orders. As suggested by Angst et ai, the energy 
difference between antiferroelectric and ferroelectric charge- 
order configuration is only ~ 3%, and it is possible to stabi- 
lize the ferroelectric configuration when the system is cooled 
in an electric field.^ The low resistivity (^ lO'^ il-m at room 
temperature) makes it very difficult to apply a static electric 
field and do the field-cooling measurements. Instead, we ap- 
plied an electric field of 20 kV/cm along the [001] direction at 
100 K, and performed scans flirough (2/3, 2/3, 3.5), (1/3, 1/3, 
0), and (0.306, 0.306, 0) peaks. The scan profiles are identi- 
cal to those shown in Fig. |2] which indicates that there is no 
observable electric-field effect with fields applied below Tco- 
At 100 K, the current is estimated to be on the order of /iA. 

To investigate the response of the system to electric cur- 
rents, we applied different currents at 360 K, cooled the sys- 
tem with the current maintained as constant, and performed 
scans at different temperatures. Results of the scans through 
the charge and magnetic peaks at 200 K with different currents 
applied are shown in Fig.|4] It is clear that for both charge and 
magnetic peaks, the peak intensities are reduced when cooled 
under electric current, and the intensity reduction increases 
with increasing current. The magnetic peak at 200 K is fully 
suppressed by a 1 mA current. To further examine the current 
effect, we plotted the charge and magnetic order peak inten- 
sity as a function of thermometer temperature with different 
currents in Fig. [Sja) and (b). 

When the sample was cooled with current, both charge and 
magnetic peak intensities were reduced. The current effects 
on the intensity clearly become more pronounced as the cur- 
rent increases. These results correlate with earlier studies 
on the electric field/current effects on transport and magnetic 



properties of LuFe204 
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and the current reducing the 



peak intensity here can be attributed to the current inducing 
breakdown of the charge order, as suggested in Refs. [TTIll2i 
However, there are also indications that this may not be as 
simple. One indication is that with 1 mA current cooling, the 
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FIG. 3: (Color Online) (a) Charge, and (b) magnetic peak intensity as 
a function of thermometer temperature, obtained under cooling with 
different currents. Inset in (b) shows temperature dependence of the 
satellite peak (0.306, 0.306, 0) intensity. Error bars represent square 
root of the total counts, and those smaller than the symbols are not 
shown. Lines through data are guides to the eyes. Arrows indicate the 
actual temperatures determined by measuring the resistance, (c) R-T 
curve measured with nearly zero (~ 1 jiA) current. The resistance 
below 170 K exceeds the maximum of the multimeter. 
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FIG. 4: (Color Online) (a) Charge (2/3, 2/3, 3.5), and (b) magnetic 
(1/3, 1/3, 0) peak measured at 200 K, after cooled with different cur- 
rents applied at 360 K. Error bars are smaller than the symbols. Lines 
through data are guides to the eyes. 



temperature reading never went below 100 K. This suggests 
that there is a significant heating power applied to the sample 
(by the current applied). Another observation is that if we re- 
move the current at 100 K, the peak intensity does not return 
to the ZFC value immediately. Instead, there is a 30 to 60 sec- 
onds lag for the (charge peak) intensity to fully recover with a 
current of 1 mA , while at the same time, the temperature read- 



1500 



1000 



^ 500 

0) 
(0 

o 

CO 

X, 



0) 

c 

0) 

~ 1000 



500 



(a) 

(2/3 2/3 3.5) 



• • ♦ • * \ 



-•-0.0 mA 

-^0.7 mA 



TV 

(1/3 1/3 0) A 
• \ 




50 



100 



150 



200 
T(K) 



250 



300 



350 



FIG. 5: (Color Online) (a) Charge, and (b) magnetic peak intensity 
as a function of calibrated temperature, obtained under cooling with 
a zero and 0.7 mA current. Error bars are smaller than the symbols. 
Lines through data are guides to the eyes. 



ing is constant. This time scale is too long for any real charge 
diffuse process to occur in these materials, and is a strong ev- 
idence that internal sample heating is playing a role. Because 
the thermometer is attached to the base of the sample mount- 
ing post, which is about 3 cm away from the sample position, 
it is plausible that there could be a large temperature gradient 
between the sample position and the thermometer location. 

In order to calibrate the sample temperature, we use the 
temperature dependence of resistance of the sample along c- 
direction as an independent measure as shown in Fig. [3jc). 
The current used to measure the resistance is small, on the or- 
der of 1 /iA, so that the condition under which the resistance 
was measured can be taken as zero-current cooling. Under 
cooling, the resistance increases continuously through Tco> 
which is consistent with the observation that below Tco- 
charge order is still short-ranged along c-axis, and disorder is 
playing an important role in the material's properties.^ Around 
250 K, the resistance drops, which is related to the magnetic 
phase transition, also suggesting a strong coupling between 
the magnetic and electrical properties in this material. Below 
250 K, the resistance increases monotonically again. 

The R-T curve provides us a good measure of the instan- 
taneous sample temperature. During the neutron diffraction 
measurements, we measured the sample's resistance by read- 
ing out the voltage across the sample while maintaining a 
constant current when cooling. Comparing the measured re- 
sistance with the R-T curve, we found that the actual sam- 
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pie temperatures under current-cooling are higher than those 
read by the thermometer The horizontal lines indicate the 
"real" sample temperature (240 K or 190 K) determined by 
the R-T curve, and the arrows indicate the temperatures read 
by the thermometer For a sample temperature of 240 K, the 
thermometer reads 240 K, 170 K, or 150 K with a current of 
0. 1 mA, 0.7 mA, or 1 mA respectively. If we correct the tem- 
perature scales of the various cooling curves using the real 
sample temperatures based on resistance readings, there is no 
field effect on the ordering. This is demonstrated in Fig.[3|a) 
and (b). On each cooling curve (with different current ap- 
plied), the data points where the real sample temperatures are 
240 K (or 190 K) are marked by arrows of different colors, 
and they indeed have the same magnetic/charge peak intensi- 
ties independent of cooling conditions. In Fig.E] we plot the 
charge and magnetic peak intensity as a function of calibrated 
temperature for cooling with a 0.7 mA current, and compare 
them with those under zero-current cooling. It is clear that 
there is no real effect if we take out the internal heating effect. 
The current-heating effect is not entirely unexpected, since 
the specific heat of LuFe204 is relatively low, ^ 0.5 J/(K g) 
at 300 K.^" With a small sample (< 1 g), the heater with a 
power on the order of 0.1 W, which corresponds to a current 
of 1 mA through the sample at 300 K, or 0.4 mA at 200 K, is 
high enough to effectively heat the sample. 

This naturally explains the unusual behaviors in LuFe204 
observed by other groups.""'^ The resistivity is low 
(^ 10'^ rt-m) in the temperature range where most measure- 
ments are carried out, so a small voltage is able to drive a 
large current through the sample and heat the sample signifi- 
cantly. The actual sample temperature in these cases will be 
higher than the readings from the thermometer The material's 
magnetization,'^ and transition temperature" '^ will then ap- 
pear to be affected by the field. The internal current heating 
can also explain the observed non-linear current-voltage be- 
havior.-'- In addition, because the resistance decreases with 



heating — if a constant voltage is applied to the sample, the 
sample is heated and the resistance lowers, which in turn in- 
creases the current further, and puts more thermal power on 
the sample, which again raises the sample temperature and 
lowers its resistance. Eventually an avalanche occurs, which 
was interpreted as the "breakdown" of the charge-orderjiisi^ 

Here we observed that the charge order in LuFe204 remains 
intact with electrical inputs — neither high electric fields ap- 
plied at low temperature, nor electric currents applied at high 
temperature can affect it. However, the charge order seems 
to be rather sensitive to magnetic field, even when no mag- 
netic order is presenti^ The fact that a charge-ordered sys- 
tem is magnetically sensitive instead of electrically sensitive 
makes LuFe204 very unusual. In other charge-ordered sys- 
tems, electric field is able to slide, or cause breakdown of the 
charge order. ^^"^'^ Apparently, from our data, it is not the case 
for LuFe204. It is likely that in LuFe204 the pinning of the 
charge order, e.g., by impurities, is stronger than that in other 
charge-ordered systems, which makes it less electrically sen- 
sitive. 

In summary, we report that the charge and magnetic order 
in LuFe204 is not affected by electric field (up to 20 kV/cm) 
or current. The observed reduction of charge- and magnetic- 
order peak intensity is due to resistive heating. Our results 
also suggest that electric field/current effects on LuFe204, as 
well as the non-linear current-voltage behavior reported else- 
where are results of internal current heating of the sample. 
This case is very similar to those observed in charge-stripe 
ordered cuprates and nickelates^^— , where an electric-field 
effect was observed on the charge order, but later entirely at- 
tributed to resistive heating effects^^. 
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